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ABSTRACT 

The thermal Sunyaev-Zel’dovich (tSZ) effect offers a means of probing the hot gas in 
and around massive galaxies and galaxy groups and clusters, which is thought to con¬ 
stitute a large fraction of the baryon content of the Universe. The Planck collaboration 
recently performed a stacking analysis of a large sample of ‘locally brightest galaxies’ 
(LBGs) and, surprisingly, inferred an approximately self-similar relation between the 
tSZ flux and halo mass. At face value, this implies that the hot gas mass fraction is 
independent of halo mass, a result which is in apparent conflict with resolved X-ray 
observations. We test the robustness of the inferred trend using synthetic tSZ maps 
generated from cosmological hydrodynamical simulations and using the same tools 
and assumptions applied in the Planck study. We show that, while the detection and 
the estimate of the ‘total’ flux (within 5rsoo) is reasonably robust, the inferred flux 
originating from within r $oo (i.e. the limiting radius to which X-ray observations typ¬ 
ically probe) is highly sensitive to the assumed pressure distribution of the gas. Using 
our most realistic simulations with AGN feedback, that reproduce a wide variety of 
X-ray and optical properties of groups and clusters, we estimate that the derived tSZ 
flux within r $oo is biased high by up to to an order of magnitude for haloes with masses 
A/500 ~ 10 13 Mq. Moreover, we show that the AGN simulations are consistent with 
the total tSZ flux-mass relation observed with Planck , whereas a self-similar model is 
ruled out. 
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1 INTRODUCTION 

Developing a detailed understanding of the formation and 
evolution of galaxies is one of the ‘holy grails’ of modern as¬ 
trophysics. It has become clear in recent years that achieving 
this goal requires developing physical models that faithfully 
capture the thermodynamic history of the gas in dark mat¬ 
ter haloes, of which only a very small fr action is able to cool 
and collapse enough t o form stars (e.g. iBaloeh et al.l[200H : 
iBudzvnski et al.ll2014l l. The problem is a challenging one, as 
the gas can be affected by a myriad of physical processes, 
such as feedback from supernovae and Active Galactic Nu¬ 
clei (hereafter AGN), radiative cooling and associated ther¬ 
mal instabilities, as well as the action of magnetic fields, 
thermal conduction, viscosity, plasma instabilities, etc. The 


situation is made even more difficult by the fact that a 
large fraction of the bary ons are ‘missing’ observationally 
(e.g. iFukueita et al .1 Pi 9981 ) . That is, only a relatively small 
fraction of the total baryonic content of the Universe (as 
inferred by measurements of the cosmic microwave back¬ 
ground, hereafter CMB, or implied by big bang nucleosyn¬ 
thesis theory combined with measurements of light elements 
at high redshift) has been mapped in any detail. Simple the¬ 
oretical calculations, as well as the predictions of detailed 
hydrodynamical simulations, suggest that a large fraction 
of the missing baryons is in the form of warm-hot (10 5 ~' 
K) gas within and surrounding dark ma t ter haloes (e.g. 
White fc Reedll978l: IWhite fc Frerikl ll Mil : I Gen fc Ostriken 
1999i : iDave et ah 200 il l, but this gas has so far been very 
difficult to detect observationally fe.g. lBregmanll2007l l. 
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X-ray observations offer one way to probe this com¬ 
ponent, but they are generally limited to relatively high- 
mass systems (clusters and X-ray-bright groups) and to only 
the central regions of nearby galaxies and low-mass groups, 
whereas the bulk of the gas/baryons are expected to be at 
much larg er ra d ius. T h e th ermal Sunyaev-Zel’dovich (tSZ) 
effect (ISunvaev fe Zeldovich|[l970l . Il972i l, which is a spectral 
distortion of the CMB background due to inverse Compton 
scattering of CMB photons off hot gas, has recently emerged 
as a new tool for probing this hot gas component. Since the 
tSZ effect depends only linearly on the gas density (as op¬ 
posed to X-ray emission, which depends on the square of 
the density) it is more sensitive to gas at large radii. Bin¬ 
ning (or ‘stacking’) large numbers of systems offers a means 
of probing the average hot gas content down to low mass. 
This technique has recently been exploited in a number 
of recent studies, such as the binning of optically-selected 


and Planck (e.g. Afshordi et al. 

2007; Bielbv & Shanks 

20071: Melin et al. 2011 

Planck Earlv Results X 201 ll: 

Planck Earlv Results XII 

201ll), as well as the binning of 


Sloan Digital Sky Survey (SDSS)-selected Luminous Red 
Galaxies using ob servations made w ith the Atacama Cos¬ 
mology Telescope (lHand et al.ll20lHl . 

Recently, the Planck collaboration reported the de¬ 
tection of significant amounts of hot gas via the tSZ ef¬ 
fect down as far as Msoo ~ 4 x 10 12 Mg (i.e. in systems 
only a factor a few more massive than the Milky Way and 
M31) by performing a binning analysis on a large sam¬ 
ple o f Locally B ri ghtest Galaxie s (LBGs) from the SDSS 
|Planck Inte rmediate Results n hereafter PIntX I; see also 
Greco et al. 1201^1 . Intriguinglv. the inferred mean relation 


between the integrated tSZ flux and halo mass is close to 
the self-similar expectation. That is, taken at face value, 
the results suggest that galaxies, groups, and clusters have 
nearly the same baryon fraction measured within rsocEl- This 
is a surprisingly result to say the least, as numerous X-ray 
studies of X-ray-bright galaxy groups and clusters have con¬ 
sistently found that groups show a marked defic it of ba ryons 
with r espect to their more ma s sive cousins (e.g. iDavid et aT 
2001 ; iGastaldello et al.l 120071 : IPratt et al.l l2009l ~ Sun et al 


20091 '). Furthermore, cosmological hydrodynamical simula¬ 
tions that invoke AGN feedback and that can reproduce 
the X-ray gas ma ss fractions of groups and clusters w i thout 
overcooling (e.g . Puchwcin ct al. 2008; Fabian et_jlJ_20^; 
McCarthy et al. 20ld : iBattaglia et al. [ l2013i: iLe Brun et alJ 
2014l ~) predict that the gas mass fractions will continue to 
steadily decline well dow n into t he individual galaxy regime. 
Therefore, if it holds, the lPInt30 result has fundamental im¬ 
plications for theories of galaxy formation and cosmic feed¬ 
back. 

However, the measurement of the tSZ signal down to 
such low masses is not a simple one by any means. Not 
only is stacking required, but strong assumptions about the 
distribution of the hot gas are required to estimate the flux 
within radii that are comparable to (or smaller than) the in¬ 
strumental beam size. Furthermore, the contribution of both 


1 rsoo is defined as the radius which encloses a mean density that 
is 500 times the critical density of the Universe at the object’s 
redshift. 


correlated and uncorrelated line-of-sight hot gas has not yet 
been well quantified. Fortunately, the impact of all these ef¬ 
fects on the recovered signal can be tested using synthetic 
observations from a suite of reasonably realistic hydrody¬ 
namical simulations. 

In the present study, we undertake a critical assess¬ 
ment of the recovery of the integrated tSZ flux (Y 500 ) as 
a function of halo mass in current observations, by using 
synthetic observations generated from a suite of state-of-the- 
art, large-volume cosmological hydrodynamical simulations 
(cosmo-OWLS). Although our methodology and tests have 
been geare d toward s a comparison with the trends reported 
recently in j T 1 1r 1 1 ■ Xfl the conclusions should be broadly appli¬ 
cable to stacked tSZ flux measurements based on matched 
filter techniques. 

The paper is organized as follows. In Section 2, we 
briefly describe the hydrodynamical simulations (Section 
2.1), our tSZ mapmaking techniques (Section 2.2), and our 
halo selection and flux recovery approach (Section 2.3). In 
Section 3, we compare the true and recovered tSZ flux-mass 
relations and quantify the difference (bias) between the two. 
In Section 4, we re-examine the pressure profiles of the hot 
gas in the simulations and design a new mass-dependent spa¬ 
tial template for the tSZ matched filter. Finally, in Section 
5, we summarize and discuss our findings. 


2 SIMULATIONS 


2.1 cosmo-OWLS 


We use the cosmo-OWLS suite of cosmolo gical hydrody¬ 
namical simulations described in detail in LeJ3run__et_alJ 
d2014j) (hereafter L14 ; see also McCart hy et alJ ~ 2014 : 
Ivan Daalen et akll2014l ; IVelliscig et ah 2014 1. They form an 
extensio n to the Overw helmingly Large Simulations project 
(OWLS: ISchave et al.llioioll . The cosmo-OWLS suite con¬ 
sists of large volume (400 h comoving Mpc on a side) 
periodic box simulations with 1024 3 dark matter and 1024 3 
baryonic particles with updated initial conditions derived ei¬ 
ther from the 7-y ear Wilkinson Microw ave Anisotropy Probe 
(WMAP ) data dKomatsu et al.l feoilh {H m , f lb, Ha, os, 


n s , h} = {0. 272, 0 . 0455, 0.728, 0.8^ 0.967, 0.704} or the 
Planck data (|Planck 2013 Results XV3 ) = {0.3175, 0.0490, 
0.6825, 0.834, 0.9624, 0.6711}. This yields dark matter and 
(initial) baryon particle masses of ~ 4.44 x 10 9 /i -1 Mg 
(« 3.75 x 10 9 h- 1 M 0 ) and ss 8.12 x 10 8 /i" 1 Mg 
7.54 x 10 8 h _1 Mg), respectively for the Planck ( WMAP7) 
cosmology. The gravitational softening is set to 4 h~ l kpc 
(in physical coordinates below 2 = 3 and in comoving co¬ 
ordinates at higher redshifts). Note that we use N ng b = 48 
neighbours for the smoothed particle hydrodynamics (SPH) 
interpolation and the minimum SPH smoothing length is set 
to one tenth of the gravitational softening. Below we present 
results based on the WMAP 7 runs only, but the results and 
conclusions are insensitive to the choice of cosmology. 

In addition to these large-volume runs, we have car¬ 
ried out additional runs with eight (two) times higher mass 
(spatial) resolution but in smaller volumes (the boxes are 
100 /i _1 comoving Mpc on a side). The motivation for this 
is two-fold: (i) to allow an exploration of the role of numer¬ 
ical convergence on the results; and (ii) to test the recovery 
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Recovery of SZ flux of dark matter haloes 3 


Simulation 

UV/X-ray background 

Cooling 

Star formation 

SN feedback 

AGN feedback 

AT hea t 

NOCOOL 

Yes 

No 

No 

No 

No 


REF 

Yes 

Yes 

Yes 

Yes 

No 


AGN 8.0 

Yes 

Yes 

Yes 

Yes 

Yes 

l 0 8.° k 

AGN 8.5 

Yes 

Yes 

Yes 

Yes 

Yes 

10 8 - 5 K 


Table 1. cosmo-OWLS runs used here and their included sub-grid physics. 


of the tSZ flux down to lower halo masses than can be re¬ 
liably probed with our larger, lower resolution box es. Note 
that the Locally Brightest Galaxy (LBG) sample of lPInt^O 
spans approximately three orders of magnitude in halo mass, 
from massive clusters with M 500 ~ 2 — 3 x 10 15 M@ down 
to individual galaxies with JW 500 ~2-3x 10 12 M®. 

The simulations were carried out with a version of the 
Lagrangian T reePM-SPH code gadget3 (last described in 
ISpringelll2005r i. which has been significantly modified to in¬ 
clude new ‘sub-grid’ physics. Starting from identical initial 
conditions (in each cosmology), the key parameters that gov¬ 
ern the most important sub-grid physics, including feedback 
from supernovae (SNe) and active galactic nuclei (AGN), are 
systematically varied. We use four of the five physical models 
described in L 14j : a non-radiative model (nocool); a model 
(ref which corresponds to the OWLS reference model) 
which in cludes prescriptions for met al-dep endent radiative 
cooling (IWiersma. Schave fc Smith! l2009aT ), stellar evolu ¬ 
tion, mass loss and chemical enrich ment ([ Wjgpnia e t al.l 
l2009bl h star formation ( Schaye fe Dalla Vecchial 20081') and 
kinetic stellar feedback ( Dalla Vecchia fe Schavel 120081 ) and 
two models (agn 8.0, which was simply called AGN in the 
original OWLS papers, and AGN 8.5) which further include 
a prescrip tion for supermassive black hole growth and AGN 
feedback dBooth fe Schavj 2009j), which is a modified version 
of the model developed bv ISpringeL Di Matteo fe Hernauistl 
(1200511 . The black holes store the feedback energy until they 
can heat neighbourin g gas particles by a p re-determined 
amount A Theat- As in lBooth fe Schavel (|2009l l. 1.5 per cent 
of the rest-mass energy of the gas which is accreted on 
to the supermassive black holes is used for the feedback. 
This results in a satisfactory matc h to the normalization 
of the b lack hole scaling relations dBooth fe Schavel 120091 : 
see also E2) which is independent of the exact value of 
A Theat- The two AGN models used here only differ by their 
value of ATheat, which is the most crucial parameter of the 
AGN feedback model in terms of the hot gas properties of 
the resulting simulated pop ulation of groups and clusters 
(IMcCarthv et al.l l201ll ; Il14 1. It is set to A Theat = 10 s K 
for agn 8.0 and ATheat = 3 x 10 s K for agn 8.5. Note that 
since the same quantity of gas is being heated in the AGN 8.5 
model as in the AGN 8.0 model, more time is needed for the 
black holes to accrete enough gas to heat the surrounding 
gas to a higher temperature. Increasing the heating tem¬ 
perature hence results into more bursty and more energetic 
feedback episodes. 

Table [T] provides a list of the runs used here and the 
sub-grid physics that they include. 

These models have b een c omp ared to a wide range o f 
ob serva tional data by both lLlj and lMcCarthv et al.l d2014lh 
InlLli we focused on the comparison to the resolved hot 
gas (e.g. X-ray luminosities and temperatures, gas fraction, 


entropy and density profiles, integrated tSZ flux) and stellar 
properties (e.g. /-band total-mass-to-light ratio, dominance 
of the brightest cluster galaxies) of local galaxy groups and 
clusters, as well as the properties of the central black hole 
and concluded that the fiducial AGN model (agn 8.0) pro¬ 
duces a realistic population of galaxy groups and clusters, 
broadly reproducing both the median trend and, for the first 
time, the scatter in physical properties over approximately 
two decades in mass (10 13 Mq < M 500 < 10 1 5 M^) and 1.5 
decad es in radius (0.05 < r/rsoo < 1-5). In iMcCarthv et ahl 
we explored the sensitivity of the thermal Sunyaev- 
Zel’dovich power spectrum to important non-gravitational 
physics and also showed that the fiducial AGN model ad¬ 
equately matches the observed pressure profiles of local 
groups and clusters (see their fig. 2). 


2.2 Theoretical thermal Sunyaev—Zel’dovich (tSZ) 
maps 


The thermal tSZ signal is characterized by the dimensionless 
Compton y parameter, defined as: 

V = [ Pedl (1) 

m,eC 2 J 

where ot is the Thomson cross-section, c the speed of light, 
m e the electron rest-mass and P e = n e fcflT e is the electron 
pressure with ks being the Boltzmann constant. The inte¬ 
gration is done along the observer’s line of sight. 

Compton y maps are generated by stacking randomly 
transformed (by a combination of translations, rotations and 
axis inversions) snapshot s along the observer’s line of sight 
(e.g. Ida Silva et al.M2000l ). The light cones extend back to 
2 = 3. Ten (sixteen) quasi-independent realisations are gen¬ 
erated for the large lower resolution (smaller higher res¬ 
olution) simulations by randomly varying the light cone 
transformations. As the methods used for the production 
of the Compton y ma ps are described in some detail in 
IMcCarthv et al.l (12014) . we will only present a brief sum¬ 
mary below. __ 

We follow the method of iRoncarelli et ahl (|2006l . 120071 ) 
and compute 


o t kBTirrii 
m e c 2 Ue,imH 


( 2 ) 


for the i th gas particle, where Ti, m; and are respectively 
the temperature, mass and mean molecular weight per free 
electron of the gas particle and m,H is the atomic mass of 
hydrogen. The contribution to the Compton y parameter by 
the i th particle is given by 

yt = Ti/ Lp ixi , (3) 

where Lp ix i is the physical area of the pixel in which the 
i th particle falls at the angular diameter distance from the 
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observer to the particle. Finally, yt is SPH-smoothed onto 
the map using the SPH smoothing kernel which was used by 
GADGET3 for the computation of the hydrodynamical forces 
and the three-dimensional smoothing length of the particle 
expressed in angular units (i.e. the smoothing length divided 
by the distance to the particle from the observer). 

We produce maps that are five degrees (1.25 degrees) 
on a side for the large lower resolution (smaller higher reso¬ 
lution) simulations. This roughly corresponds to the angular 
size of the 400 h comoving Mpc (100 h comoving Mpc) 
box at z = 3.0. The maps have an angular pixel size of 2.5 
arcseconds. 

As discussed in lMcCarthv et ahl (I2014I ). because the re¬ 
alisations are produced using the same simulations, they are 
only quasi-independent. At high redshift, the light cone sam¬ 
ples most of the volume of simulation and thus the different 
realisations contain many of the same structures (though 
at randomised locations). At lower redshift, however, the 
cones sample only a small fraction of the simulated vol¬ 
ume and the different realisations are therefore effectively 
independent. This is relevant for the present study, which 
is focused on the recovery of the tSZ flux of relatively local 
(z ^ 0.4) dark matter haloes. The rationale behind using 
several quasi-independent realisations is to mitigate the im¬ 
pact of cosmic variance by significantly increasing the sam¬ 
ple size. 

2.3 Halo selection and tSZ flux recovery 

The tSZ signal of individual galaxies, groups, and low-mass 
clusters cannot normally be detected with Planck. Large 
numbers of systems must therefore be binned (‘stacked’) 
to measure the mea n re lation between haloes and tSZ flux 
down to low masses. IPIntXj undertook such a binning anal¬ 
ysis of some ~ 260 000 locally brightest galaxies selected 
from the SDSS. They defined their locally brightest galaxy 
(hereafter LBG) sample as all the galaxies with 2 > 0.03 and 
r < 17.7 which are brighter than any other galaxy within a 
projected distance of 1.0 Mpc and with a redshift difference 
smaller than 1,000 km s -1 . The sample was derived from 
the spectroscopic New York LTniversity Value Added Galaxy 
Catalogue (NYU-VAGCfl which is based on the seventh 
data release of the Sloan Digital Sky Survey (SDSS/DR7; 
lAbazaiian et al.ll2009l ') and covers 7966 square degrees. The 
NYU-VAGC also provides, amon g other properties, stellar 
masses which were computed by IBlanton fe Roweia (|2007l l 
by fitting stellar pop ul ations to t he five-band SDSS photom¬ 
etry assu ming a lChabriei] (|2003l l initial mass function. 

IPIntXj derived the mean tSZ flux in bins of stellar mass. 
To deduce the more fundamental relation between tSZ flux 
and halo mas s, th e y use d the semi-analytic galaxy forma¬ 
tion model of I Guo et al.l (1201 il l . which was tuned to closely 
match the observed luminosity and stellar mass fu nct ions o f 
SDSS galaxies in a WMAP7 cosmology (Guo et ahll2013l) . 
When deriving the relation between tSZ flux and halo mass, 
the effects of scatter in the stellar mass-halo mass relation, 
as well as (minor) contamination of the LBG sample by 
satellite systems, are accounted for using the semi-analytic 
model. 

2 http://sdss.physics.nyu.edu/vagc/ 


Our aim is to test the accuracy of the recovered tSZ 
flux-halo mass relation by analysing synthetic maps in a 
way that is faithful to that done for the real data. We place 
our emphasis on the recovery of the tSZ flux, rather than 
testing the conversion of the tSZ flux-stellar mass relation 
into an tSZ flux- halo mass relation. The latter will depend 
on the models correctly capturing the intrinsic scatter in the 
stellar mass-halo mass relation (i.e. in any given stellar mass 
bin, there will be a range of halo masses and since the SZ flux 
has a steeper than linear scaling with halo mass, the mean 
signal will be driven by the higher halo mass systems in the 
bin). The scatter can be constrained to some degree by re¬ 
quiring the models to match, for example, the galaxy stellar 
mass function and/or the clustering of gala xies in b ins o f 
stellar mass (as is the case for the model of lGuo et al.ll201ll . 
used in the Planck LBG study), but it is probably fair to 
say that the intrinsic scatter issue is still under considera¬ 
tion. This uncertainty in the stellar mass-halo mass relation 
could therefore potentially further comp licate the physical 
interpretation of the results presented in IPIntXj . 

We first construct halo catalogues corresponding to the 
maps described in Section m\ using a standard friends-of- 
friends (FoF) algorithm run on the snapshot data. The cat¬ 
alogues contain the positions on the map, the angular size 
#500 (= rsoo/dA, where dA is the angular diameter distance 
to the halo) and the halo mass M 500 of all the haloes with 
z < 0.4 and AL 500 > 2 x 10 13 Mq (A /500 > 10 12 Mq) for 
the large lower resolution (smaller higher resolution) sim¬ 
ulations. The redshift and mass thresholds h ave been cho¬ 
sen to roughly reproduce the bounds of the IPIntXll LBG 
sample. To produce ‘dirty’ maps with similar characteristics 
to the real data, the raw maps are downgraded from their 
original 2.5 arcseconds resolution to a 0.83 arcminute reso¬ 
lution and the effects of primary CMB, of the Planck beams 
and their associated noise are added to obtain six syn¬ 
thetic observations at the frequencies of the Planck HFI 
instrument (100, 143, 217, 353, 545 and 857 GHz). Note 
that in producing the dirty maps, we have neglected dust 
and radio point source contamination, which can be signifi¬ 
cant for real observations. However, this omission is at least 
partially compensated for by the fact that the overlapping 
SDSS -Planck survey area (7966 square degrees) greatly ex¬ 
ceeds the survey area that can be simulated with current 
self-consistent cosmological hydro simulations — the cosmo- 
OWLS maps, which are some of the largest ever produced, 
are still only 25 square degrees each. In practice, we produce 
as many (primary CMB+instrumental) noise realisations of 
the simulation maps as are needed to obtain statistical er¬ 
ror bars on the derived tSZ flux-hal o mass r elation that 
are comparable to the ones reported in|PIntX|, namely two 
(twenty-five) for the large lower resolution (smaller higher 
resolution) simulations. 

Following the production of the dirty ma ps, we the n ap¬ 
ply t h e same multi - frequ ency matched filter (iHerranz et al.l 
120021 ; iMelin et al.l I2006I . MMF) that was used by the 
Planck collaboration. The MMF yields an estimate of the 
tSZ flux within 5rsoo, which is then converted using a con¬ 
stant conversion factor into the tSZ flux within rs 00 (see 
discussion below), as characterized by the value of its spher¬ 
ically integrated Compton parameter within r^oo'- 

d A {z) 2 Y 500 = f p e dV (4) 

m e c 2 ./ 
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where cIa(z) is the angular diameter distance and the inte¬ 
gration is done over a sphere of radius rsoo. Throughout the 
paper, unless otherwise stated, we use the quantity 

= (5 ) 

where E(z) = H(z)/Ho = + z) 3 + S!a gives the 

redshift evolution of the Hubble parameter H(z) in a flat 
ACDM Universe. I 500 corresponds to the intrinsic tSZ flux 
self-similarly scaled to z = 0 and scaled to a fixed angular 
diameter distance. 

The MMF is optimised in both frequency and angu¬ 
lar spactjf] by assuming the known frequency dependence of 
the thermal tSZ effe ct and th e ‘un iversal pressure profile’ 
llArnaud et alj l20ld . hereafter I Aldl ') , derived from a com¬ 
bination of X-ray obser vations of the XMM-N ewton REX- 
CESS cluster sample llBohringer et al.l 1200711 at radii of 
r < r 500 and hydrodynam ical sim ulations at radii of rsoo < 
r < 5r5oo ■ Consistent with lPIntXl the MMF is run in a non¬ 
blind mode using the positions and sizes, # 500 , from the halo 
catalogue. The MMF then gives a measure of the strength of 
the tSZ signal lsoo(*) and its associated measurement un¬ 
certainty ag 500 (i) for the halo surrounding the i th galaxy. 
The measurement uncertainty takes into account the statis¬ 
tical uncertainties due to astrophysical (e.g. primary CMB) 
and instrumental noise, but not the uncertainties due to 
halo modelling (e.g. shape of the pressure profile, size). Note 
that the tSZ spectral function was not integrated over the 
Planck bandpasses for each frequency for both constructing 
the ‘dirty’ maps and creating the matched filter for each of 
the Planck frequencies. This should have no impact upon 
the results as it was done at both stages. 

It is important to note here that the LBGs are at best 
only marginally resolved by Planck. Therefore, in practice, 
the tSZ flux is actually measured within the larger aperture 
of 5#5oo; 

050500 

Y sr 500 = / 2t v6y(d)de (6) 
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Figure 1. Comparison of the Y 50 Q—M 500 relation inferred by 
IPIntXll with the intrinsic (true) z = 0 relation from the cosmo- 
OWLS runs. The fille d cyan c ircles with error bars represent the 
observational data of IPIntxfl . The solid and dashed curves (red, 
orange, blue and green) represent the mean and median inte¬ 
grated tSZ flux-Msoo relations in bins of M 500 f or the d ifferent 
simulations, respectively. Taken at face value, the IPIntXj results 
favour a close to self-similar Y 5 oo - ^500 relation. However, Y 500 
is not directly measured and synthetic observations are therefore 
required for proper comparison to the inferred trend (see Fig. [ 3 . 


N b 


a b 2 = X 1 / 5 ioo(*)’ 


(8) 


where Nb is the number of galaxies in bin b. The integrated 
tSZ flux is often normalized by the self-similar i ntegr ated 
tSZ flux-mass relation, as given in appendix B of Aid : 


The flux within the spherical aperture rsoo is then com¬ 
puted assuming the spatial template used in the MMF (the 
universal pressure profile), leading to a conversion factor 
l 5 00 = y'5 I . 500 /1.796. (It is assumed that there is no flux orig¬ 
inating from beyond 5rsoo and therefore Y 5 g 500 = lsr 500 ■) It 
is perhaps an obvious statement that if in reality the adopted 
spatial template poorly describes the pressure distribution 
of the hot gas, this will lead to a systematic error in the re¬ 
covered flux within rsoo- We will return to this point below. 

Finally, for both simulations and observations, the tSZ 
flux is binned by mass with the bin -average flux and the cor¬ 
responding uncertainty g i ven by dPIanck Early Results Xt 
IPIanck Early Results xfl IPIntXj) 


(^500)5 


E ^noo (»)/<„„(») 
ESi V^Sboo (*) 


and 


(7) 


3 It maximises the signal-to-noise ratio of objects which follow 
the assumed spectral and spatial templates. 


^50o,aio — 9.07 x 10 


M 500 


.3 x 10 14 /itoM© 


5/3 


h 7 f arcmin 2 (9) 


3 THE NEED FOR SYNTHETIC TSZ 

OBSERVATIONS 

We begin by comparing, in Fig.p] the lPIn txi results to the 
true spherically integrated tSZ flux-mass (Y 500 — A/ 500 ) rela¬ 
tion of the simulated galaxy, group and cluster populations 
for the four physical models at z = 0. The solid and dashed 
curves respectively represent the mean and median rela¬ 
tions in bins of M 500 for the different simulations. The high- 
resolution simulations are used below M 500 = 2 x 10 13 Mg 
and the production runs above that threshold (see Appendix 
A). The filled cyan circles with error bars represent the ob¬ 
servational data of lPIntXl 

Focusing first on the simulations, there is a visible steep¬ 
ening of the I)ioo — A /500 relation in the galaxy and group 
regime (M 500 < 10 14 Mg) with the inclusion of AGN feed¬ 
back. This is due to the fact that the AGN lower the gas 
fractions, and hence the integrated Compton y values, of 
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Figure 2. Comparison of the Y^oq-M^qq relation inferred by IPIntXj with the recovered relation from synthetic tS Z ob serv ations of 
the cosmo-OWLS runs. The blue dashed and cyan dot-dashed lines correspond to the best-fitting scaling relations of I Aldl and PlntXI, 
respectively. The solid and dotted black curves correspond to the unweighted and weighted mean true tSZ flux M 500 relations in bins 
of M 500 for the diffe rent sim ulations, respectively. The filled cyan circles and empty black diamonds with error bars correspond to the 
observational data of PlntXI and the results of the synthetic obser vations of the simulated maps, respectively. The integrated tSZ fluxes 
have been normalized using the best-fitting scaling relation of Aid . The different panels correspond to the different physical models. The 
MMF recovers the intrinsic Y 500 AAoo well for models that neglect efficient feedback (NOCOOL and ref). However, a significant bias is 
present in the recovered relation for models which include AGN feedback. The bias increases strongly with decreasing halo mass and 
with increasing AGN heating temperature (i.e. more violent feedback events). 


low-m ass haloes fe.g. IPuchwein et al.ll2008l : iMcCarthv et al.l 
l20ld ; L14) by ejecting gas from their high-redshift progeni¬ 
tors llMcCarthv et all 201 ll ). 

Interes tingly, it is apparent that, taken at face value, 
the IPIntXll results appear to favour an approximately self¬ 
similar Y 500 -M 500 relation, as obtained in non-radiative sim¬ 
ulations (such as NOCOOL; red lines). This is consistent with 


the conclusions of IPIntXll As already discussed in Section 
1 , this is a surprising result, given that numerous previ¬ 
ous studies have demonstrated that the self-similar model 
is s tron gly ru led ou t by X-ray observations. In L14 (see 
also IMcCarthv et al.ll20ld) , for example, we showed that the 
fiducial AGN model (agn 8.0) reproduces a wide variety of 
X-ray and optical measurements of local groups and clus- 
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ters (including their gas fractions) while the NOCOOL model 4 ROLE OF DEVIATIONS FROM THE 
is strongly disfavoured. UNIVERSAL PRESSURE PROFILE 


As discussed in Section however, the quantity Vsoo 
is not directly measured with Planck. Instead, what is mea¬ 
sured is an estimate of the integrated flux within a 2D aper¬ 
ture 50500, and a scaling factor which depends on the as¬ 
sumed distribution of the hot gas must be applied to esti¬ 
mate the tSZ flux within the spherical radius rsoo- Further¬ 
more, hot gas along the line-of-sight can potentially bias the 
2D flux. The magnitude of these effects can be quantified us¬ 
ing synthetic observations of the type described in Section 
2.3. 


With this in mind, in Fig. 0 we compare the recovered 
tSZ flux—A /500 relation from our synthetic (‘dirty’) maps 
(empty black diamonds) with the mean true tSZ flux-A/500 
relation, both unweighted (solid lines) and weighted (dotted 
lines; as given b y equati on 0 ) , as well as with the observa¬ 
tional data from lPIntXIl (filled cyan circles with error bars). 
To more clearly emphasize the level of agreement (or lack 
thereof) between the recovered and true relations and the 
observational data, we n orma lize the tSZ flux using the best¬ 
fitting scaling relation of lAldl (see equation [9]). The different 
panels correspond to the different physical models: from NO¬ 
COOL (top left) to AGN 8.5 (bottom right) through ref (top 
right) and AGN 8.0 (bottom left). 


Focusing first on the top two panels of Fig. [2] corre¬ 
sponding to the NOCOOL and REF simulations, we find that 
the MMF recovers the true relations relatively well, as de¬ 
duced by comparing the empty black diamonds with the 
black curves. (The only exception to this is the lowest-mass 
bin, where the recovered flux is significantly higher than 
the intrinsic one, which we discuss below in more detail.) 
It is interesting to note that the recovery is best for the 
NOCOOL model, which has a nearly self-similar behaviour, 
in the sense that the shape of the pressure distribution of 
the hot gas is nearly independent of halo mass. It is worth 
re-iterating that the (empirical) universal pressure profile 
of A10 was derived from a sample of very massive clusters, 
where the high mass basically ensures that non-gravitational 
physics (e.g. AGN) has a relatively minor effect on the tSZ 
flux. Thus, the fact that the MMF recovers the true relation 
well for a simulation where the gas approximately follows 
the universal pressure profile on all mass scales may not be 
that surprising. In fact, it is reassuring. 


Moving on to the bottom two panels of Fig. 0 cor¬ 
responding to the two AGN models, we find a significant 
offset between the recovered and true Y 500 -M 500 relations. 
The amplitude of the bias increases strongly with decreas¬ 
ing halo mass and with increasing AGN heating tempera¬ 
ture. Notably, the bias reaches nearly an order of magnitude 
at the lowest halo masses. The fact that the models with 
AGN much more closely reproduce the X-ray and optical 
properties of real groups and clusters than the NOCOOL and 
ref simulations provides a strong impetus to understand 
the origin of the bias and to offer a means for correcting for 
it. 


As the magnitude of the bias between the recovered and 
true V 500 -A /500 relations in Fig.0cvidently depends strongly 
the nature of non-gravitational (sub-grid) physics included 
in the simulations, this suggests its origin is tied to devia¬ 
tions of the hot gas distribution in the simulations from that 
adopted in the matched filter. We therefore explore this pos¬ 
sibility in detail below. An alternative source of bias at the 
low-mass end is from source confusion (i.e. hot gas along 
the line of sight). We explore the role of source confusion 
in Appendix B, demonstrating that, while it is expected to 
significantly increase the scatter in the recovered fluxes of in¬ 
dividual haloes, source confusion does not significantly bias 
the recovered mean tSZ flux-halo mass relation. 

The ‘universal pressure profile’ of A10 is described 
well by a so-called general ized NFW (GNFW) model 
and was first applied by iNagai et al.l (2007) to de¬ 
scribe the pressure distribution of the hot gas in clus¬ 
ters. It was subsequently shown to provide a rel¬ 
atively good description to large samples of (high- 


mass) 

X-ray groups and clusters (e.g. Mroczkowski et al. 

20091: 

2011; 

Arnaud et al. 201Q 1 : Plagge et al. 201$” Sun et al. 

Planck Intermediate Results Vi: Savers et al. 2013: 

McDonald et al. 2014). The model has five free parameters: 


P(r) = _ Po _ 

R 500 (c5oor/r 50 o) 7 [l + (c5oor/r5oo)“P -7)/a ’ 

where 


( 10 ) 


P500 = n ej 5 oofcflT 5 oo ( 11 ) 

500/;, Pcrit (^) l-tm-p G A /500 (12) 

/.i e mjr2r5oo 

fb = Plb/Plm is the universal baryon fraction, /r and are 
the mean molecular weight and the mean molecular weight 
per free electron, respectively. 

A10 fitted the generalised NFW profile to the REX- 
CESS sample of X-ray clusters to constrain four of the 
five p arameters (Po, C500, 7, and a). Following iNagai et al.l 
{ 2 OO 3 ), A10 fixed the value of j3 (the external slope) to 5.49, 
the value which best describe s the large-radii behaviour of 
the simulations of lNagai et al.l (120071 4 . We note here that the 
pressure distribution of the gas as inferred by X-ray obser¬ 
vations is only constrained out to r ~ rsoo • Hydrodynamical 
simulations have therefore been relied on by previous studies 
to constrain the outer slope. However, recently resolved tSZ 
observations of local clusters with Planck have now allowed 
for direct measurements of the pressure profile out to several 
rsoo and indicate that the external slope may be somewhat 
fla tter (/3 ~ 4) than predicted by some previous simulations 
(jpianck Intermediate Results Vl l. Interestingly, our cosmo- 
OWLS models predict outer pressure d istributions th at are 
in very good agreement with lPlanck Intermediate Results Vl 
(see Fig. 0below, comparing the simulated profiles with the 
universal pressure profile assuming j3 = 5.49). One possible 
reason for the difference in pressure profiles between cosmo- 
OWLS and the simulations used by A10 is that the previous 
simulations lacked feedback from AGN, which is crucial for 
preventin g overcoolin g, which should lead to reduced pres¬ 
sures. iBattaglia et a l.l (i2f)12h . who use the AGN simulations 
of lBattaglia et al.l ( 2010h . have come to similar conclusions. 
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AGN 8.0 

1.001 . ’ ”1 



0.1 1.0 0.1 1.0 

r / r 500 


Figure 3. Mass-weighted pressure profiles at z = 0 for the AGN 8.0 model. The eight panels correspond to eight different mass bins 
(increasing from top left to bottom right). The bottom part of each panel shows the residuals of the best-fitting profiles. The black 
filled circles with error bars correspond to the median profile in the corresponding mass bin and the error bars encompass 68 per cent 
of the systems. Th e soli d blue, dashed blue, solid cyan, orange, magenta and r ed lin es correspond to the universal pressure profile of 
the appendix B of I Aid , the empirical universal pressure profile of section 5 of I Aid , the best-fitting GNFW functional form with all 
the parameters left free to vary, the best-fitting GNFW form but with the concentration now a power-law of mass, the best-fitting 
GNFW form but with the normalization now a power-law of mass, and the best-fitting GNFW form but with the concentration and 
normalization now power-laws of mass, respectively. The shape of the pressure profiles is quite strongly mass-dependent and two of the 
GNFW coefficients must be made mass-dependent in order to obtain a reasonable fit (solid red curve) over the whole radial and mass 

ranges when AGN feedback is included. __, 
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Simulation 

Median or mass-weighted 

Po,o 

a 

P 

7 

C500,0 

5 

e 

REF 

Mass-weighted 

0.528 

2.208 

3.632 

1.486 

1.192 

0.051 

0.210 

REF 

Median 

0.694 

1.489 

4.512 

1.174 

0.986 

0.072 

0.245 

AGN 8.0 

Mass-weighted 

0.581 

2.017 

3.835 

1.076 

1.035 

0.273 

0.819 

AGN 8.0 

Median 

0.791 

1.517 

4.625 

0.814 

0.892 

0.263 

0.805 

AGN 8.5 

Mass-weighted 

0.214 

1.868 

4.117 

1.063 

0.682 

0.245 

0.839 

AGN 8.5 

Median 

0.235 

1.572 

4.850 

0.920 

0.597 

0.246 

0.864 


Table 2. Results of the pressure profile fitting. 


We note that iNagai et al.l (j2007t ) justified the choice 
of this particular functional form by the fact that the 
gas pressure distribution is primarily determined by the 
gravitationally dominant dark matter, in which the gas 
is in hydrostatic equilibrium, and whose density has been 
sh own to approximately fol low the NFW profile (e.g. 
[Navarro, Frenk fc Whitefl997l ) . While gas tracing dark mat¬ 
ter may be an acceptable approximation for relatively mas¬ 
sive systems where the effects of feedback from AGN and 
supernovae are more minor, X-ray observations show the 
gas is signi ficantly more extende d (less cuspy) in low-mass 
groups fe.g. lVikhlinin et all2006h . This behaviour is evident 
in our simulations with AGN feedback as well. At the very 
least, therefore, we anticipate that some of the parameters 
of the GNFW profile will need to be functions of halo mass, 
rather than fixed constants (as assumed in A10), in order 
to describe the pressure distribution of the hot gas over the 
full range of masses in the cosmo-OWLS AGN models. 


With the above in mind, we show, in Fig. [3] the di¬ 
mensionless mass-weighted pressure profiles for the AGN 8.0 
model in eight different mass bins (increasing from top left 
to bottom right; the same bins were used for binning the 
tSZ flux in Fig. 0. We normalize the radii by rsoo and the 
pressures by P 500 • Note that in order to further reduce the 
dynamic range of the y-axis, we plot P/P 500 (r/rsoo) 2 - Note 
also that the X-ray pressure profiles of the simulated groups 
and clusters from cosmo-OWLS have alr eady been compared 
to the X-ray observa tions of lSun e t al. d 2011 ) (groups) an d 
lArnaud et al.l (120101 ) ( clusters) bv iMcCarthv et al.l d2014h . 
who found that the fiducial AGN model (agn 8.0) under 
consideration here reproduces the pressure profiles of both 
groups and clusters well. 


In Fig. EO the black filled circles with error bars corre¬ 
spond to the median profile in the corresponding mass bin 
and the error bars encompass 68 per cent of the systems. 
The solid blue, cyan, orange, magenta and red lines corre¬ 
spond to the universal pressure profile of the appendix B 
of IaToI. the best-fitting GNFW functional form with all five 
parameters left free to vary, the best-fitting GNFW form 
but with the concentration now a power-law of mass of the 
form C500 = C 5 oo,o(Af 5 oo/ 10 14 M 0 ) s , the best-fitting GNFW 
form but with the normalization now a power-law of mass 
of the form Po = -Po,o(Af 5 oo/ 10 14 M@) e , and the best-fitting 
GNFW form but with the concentration and normalization 
now power-laws of mass, respectively. The bottom part of 
each panel shows the residuals of the best-fitting profiles. 
The functional forms are fitted simultaneously to the eight 
median profiles with their error bars (i.e. to the filled black 
circles with error bars). The best-fitting parameters for the 
final model which has both the concentration and normal- 


AGN 8.0 



Figure 4. Comparison of the Vsoo - -^500 relations when a new 
mass-dependent spatial template based on the AGN 8.0 model is 
used in the matched filter. The blue dashed and cyan d ot-da shed 
lines co rrespond to the best-fitting scaling relations of [Aia and 
IPIntXj , respectively. The solid and dotted black curves corre¬ 
spond to the unweighted and weighted mean true tSZ flux-Msoo 
relations for the AGN 8.0 model, respectively. The filled cyan cir¬ 
cles and empty black di amonds with error bars correspond to 
the observational data of lPIn tx] and the results of the synthetic 
observations (using the new template) of the simulated maps, 
respectively. The integrated tSZ fluxes h ave been normalized us¬ 
ing the best-fitting scaling relation of lAldl . The recovered and 
true Y 500 -M 500 approximately agree, implying the majority of 
the bias present in Fig. [3] is due to the use of an inappropriate 
(fixed) spatial template. 


ization varying as power-laws of mass for the four cosmo- 
OWLS models used here and for both mass-weighted and 
median pressure profiles are listed in Table [2] We note that 
some of the parameters are strongly degenerate with each 
other, but this is unimportant for our purposes since all that 
we require is a parametric model that reproduces simulated 
pressure profiles well. 

From Fig. [3] it is clear that the shape of the pressure 
profiles is indeed quite strongly mass-dependent (compare 
the black filled circles to the solid blue curves, which rep¬ 
resent the fixed-shape universal pressure profile) and it is 
necessary to make two of the GNFW coefficients, the nor¬ 
malization and the concentration, mass-dependent in order 
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to obtain a reasonable fit (solid red curve) over the whole 
radial and halo mass range when AGN feedback is included. 
Note that iBattaelia et al.l 1 2012 ) similarly found that the 
normalization and concentration as well as the outer slope 
/3 needed to be not only mass-dependent but also redshift- 
dependent in orde r to pr ovide a decent fit to the simulations 
of iBattaglia et al . I <l2010l) which include AGN feedback. 

We now proceed to alter the template adopted in the 
MMF, replacing the universal pressure profile with the tem¬ 
plate derived above (the solid red curve in Fig. [3j and we 
repeat the recovery of tSZ fluxes from the synthetic maps. 
In Fig. H we compare the results of the synthetic observa¬ 
tions of the simulated maps for z < 0.4 based on the new 
template (empty black diamonds) with the mean spherically 
integrated (true) tSZ flux-Msoo relation (black curves), as 
well as the observational data (filled cyan circles with error 
bars). Note that when using the best-fitting pressure pro¬ 
file taken from the simulation as a spatial template in the 
matched filter, we numerically compute the conversion fac¬ 
tor from Y$ r500 into I 500 using the new template profile for 
consistency. In Appendix C, we show the derived conversion 
factor as a function of halo mass using the new template and 
compare it with the true conversion factor derived directly 
from the simulations. 

Quite remarkably, the bias between the recovered and 
true values of Y 500 for the AGN 8.0 model in Fig. [2]is largely 
removed. Hence, the majority of the bias between the recov¬ 
ered and true tSZ flux within r 500 is due to the use of an 
inappropriate spatial template in the matched filter. 

Note that even with the new template, some residual 
bias remains, implying the template does not fully capture 
the detailed structure of the hot gas (e.g. halo-to-halo scat¬ 
ter may be significant, gas clumping, asphericity, etc.). In¬ 
deed, as we show in Appendix C, the conversion factor (from 
Isrsoo to V 500 ) trend derived from the best-fitting mass- 
dependent spatial template does not perfectly match the 
true conversion factor relation from the simulations. One 
may be able to derive a spatial template that better repro¬ 
duces the true conversion factor from the simulations, but 
this would likely come at the expense of the template provid¬ 
ing a poorer fit to the spherically-symmetric pressure profile. 
Furthermore, even perfect knowledge of the conversion fac¬ 
tor does not guarantee a bias-free estimate of I 500 , as the 
estimate of ‘total’ flux (Ysrsoo) may also be biased at some 
level, as we discuss immediately below. 

It is of interest to determine whether the removal of the 
bias (in Fig. II with the use of the new template is due en¬ 
tirely to the improved conversion in the flux measured within 
5r5oo to that within rs 00 , or if the ‘total’ flux within 5rsoo 
itself has also been affected. That is, is the flux measured 
within the larger radius 5r5oo unbiased? In the left panel 
of Fig. [5] we show the unweighted mean true Ys r500 -Msoo 
relation for the AGN 8.0 simulation (solid lines), the observa¬ 
tional data (filled cyan circles with error bars), the recovered 
relation from synthetic observations of the simulated maps 
for « < 0.4 when using the new halo mass-dependent spatial 
template (empty black diamonds), and the recov ered rela¬ 
tion when using universal pressure profile oflAlol (standard 
MMF; filled blue circles), respectively. Even when measured 
within 5r5oo, the recovered fluxes are still biased at a rela¬ 
tively low level (~ 10 — 20 per cent) with respect to the true 
Ysr 500 — Msoo (solid line) for the AGN models. Using our new 


mass-dependent template improves the situation somewhat, 
but does not eliminate the bias altogether. 

For the sake of consistency and completeness, in the 
right panel of Fig. [5] we compare the mean Y 5 ,. 500 — M 500 
relation for the synthetic observations using the standard 
MMF on the simulated maps for the different physical mod- 
els (c oloured solid curves) with the observational data of 
IPIntXj (filled cyan circles). This represents a like-with-like 
comparison between all the models and the observational 
data, in the sense that we are directly comparing what has 
actually been measured and the measurements for both the 
observations and the simulations have been made in a consis¬ 
tent manner (even if the standard MMF template has been 
shown to poorly reproduce the hot gas in some of these mod¬ 
els). We conclude from this comparison that t here is v ery 
good agreement between the trend measured by lPIntxH and 
that predicted by the fiducial AGN (agn 8.0) model which 
reproduc es th e X-r ay properties of local g roups and clusters 
well (see lL14l and iMcCarthv et alJ 12014 1. In contrast, the 
self-similar (nocool) model is excluded on the basis that 
it overpredicts the observed tSZ flux at low to intermediate 
halo masses. 


5 SUMMARY AND DISCUSSION 

We have used synthetic observations produced from the 
cosmo-OWLS suite of cosmological hydrodynamical simula¬ 
tions to: (i) perform a consistent (like-with-like) comparison 
with recent observations; and (ii) to test the reliability of 
the observationally-inferred integrated tSZ flux-halo mass 
relation. While our methodology and tests have been geared 
towards a comparison wit h the trends recentl y rep orted in 
IPlanck Intermediate Results XI (hereafter [PlntXIj), the re¬ 
sults and conclusions should be broadly applicable to tSZ 
flux measurements based on matched filter techniques. 

From the analysis presented here, we reach the following 
conclusions: 

(i) Taken at face value, the results of PIntXI favour a 
close to self-similar Y 500 -M 500 relation (Fig. UJ. Assuming 
that the gas has a temperature that is close to the halo virial 
temperature, such a trend would imply that the haloes of 
galaxies, groups and clusters all have a gas fraction within 
rsoo that is close to the universal baryon fraction and inde¬ 
pendent of mass. However, X- ray observations of local X- 
ray-bright galaxy groups (e.g. I Sun et ahl 120091 ) show that 
they are missing a large fraction of their baryons, indicat¬ 
ing some inconsistenc y betw een the tSZ- and X-ray-derived 
results, as noted bv IPIntXll 

(ii) We find that the multi-frequency matched filter 
(MMF) used by the Planck collaboration recovers a nearly 
unbiased Y 500 -M 500 relation for models which neglect effi¬ 
cient feedback (NOCOOL and ref ; see top panels of Fig. [5J. 
This is perhaps not too surprising (it is reassuring), since for 
these models the hot gas distribution has a shape similar to 
the universal pressure profile that is nearly independent of 
halo mass. 

(iii) However, we find a significant bias in the recovered 
Y 500 -A /500 relation (see bottom panels of Fig. HI) for models 
which invoke AGN feedback, which is necessary to reproduce 
the X-ray and optical properties of local groups and clusters 
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Figure 5. Compa rison of t he Yr >,mn Myxi relations. Left: The blue dashed and cyan dot-dashed lines correspond to the best-fitting 
scaling relations of I Aldl and IPIntXj . respectively. The solid black curve correspond to the unweighted mean true l" 5 r 500 —M 500 relation 
for the AGN 8.0 simu lation . The filled cyan circles, empty black diamonds and filled blue circles with error bars correspond to the 
observational data of [Plntxl and the results of the synthetic observations of the simulated maps with the extraction made using the 
new mass-dependent spatial template (derived from the AGN 8.0 run) and the standard template of I Aldi respectively. Overall, the level 
of bias present in the recovered flux within oryio is relatively modest (10 — 20 per ce nt), inde pendent of the choice of spatial template. 
Right: A like-with-like comparison of the cosmo-OWLS runs with measurements of IPIntxl . The solid curves (red, orange, blue and 
green) represent the mean relations for the different physical models using the synthetic observations with the standard template. When 
comparing observed quantities (i.e. Y 5 r500 instead of Y 500 ), the Planck observations actually exclude self-similar behaviour (nocool) 
rather than support it. The fiducial AGN model (AGN 8.0), which reproduces the X-ray and optical properties of local groups and clusters, 
is consistent with the measured tSZ flux—mass trend. 


(see iMcCarthv et al.l l20ld : iLe Brun et al.l l20l4) . The am¬ 
plitude of the bias increases strongly with decreasing halo 
mass, reaching nearly an order of magnitude overestimate 
in F 500 at halo masses below ~ 10 13 Mg. 

(iv) We have shown that the vast majority of the bias 
originates from the assumption of a fixed spatial template 
(the universal pressure profile), which becomes an increas¬ 
ingly poor description of the hot gas at low masses in our 
models with AGN feedback (Fig. El- 

(v) When a mass-dependent spatial template that de¬ 
scribes the pressure profiles in the AGN models well is used 
(Fig. E]and Table EJ, the fluxes recovered for that simula¬ 
tion are nearly unbiased with respect to the true 1500 — M 500 
relation of that simulation (Fig. El ■ 

(vi) Encouragingly, fluxes measured within the larger 
aperture 5rsoo (i.e. F 5 r 500 ) only show modest levels of bias 
(~ 10 per cent), independent of the detailed shape of spa¬ 
tial template (left panel of Fig. El or the nature of the feed¬ 
back implemented in the simulations (i.e. whether or not 
AGN feedback is included). Thus, the reported detection of 
a large fraction of the ‘missing baryons’ in the form of hot 
halo gas by the PIntXI appears to be quite robust, even if 
th e gas is more e xtended than previously thought (see also 
iGreco et al.ll2014l ). 

(vii) A consistent like-with-like comparison between the 
tSZ fluxes measured within the larger aperture 5rsoo demon¬ 
strates that the fiducial AGN model is in excellent agree¬ 


ment with the trends measured in |PIntx3 . whereas the self¬ 
similar (nocool) model actually overpredicts the tSZ flux 
at low to intermediate halo masses (right panel of Fig. EJ- 

The analysis of our AGN models indicate that the es¬ 
timated tSZ fluxes within rsoo may be significantly biased 
(high) at low to intermediate masses. However, it should be 
borne in mind that the estimate of the bias is quite sensi¬ 
tive to the nature of the feedback implemented in the sim¬ 
ulations. We have used our current ‘best guess’, as judged 
by previous comparisons of these models with resolved X- 
ray, tSZ, and optical observations (L14). Independent tests, 
however, are certainly well worthwhile. In particular, an ob¬ 
vious prediction of our results is that future high-resolution 
tSZ observations (e.g. SPTpol, ACTpol, SPT-3G, AdvACT, 
NIKA), which are capable of directly measuring the tSZ flux 
within 05 oo, ought to see a steepening in the Y 500 --M 500 re¬ 
lation below M 500 ~ 10 14 Mq due to gas ejection by AGN. 

An alternative test is to stack X-ray observations of a 
large optically-selected sample (such as the SDSS LBG sam¬ 
ple in PI ntXI) and to directly measure the hot gas proper¬ 
ties within rsoo in b ins o f mass f or co mparison to the tSZ- 
derived results. lAnderson et al.l (120151 ) have recently under¬ 
taken such an exercise using the ROSAT All-Sky Survey 
(RASS). They derive the stacked X-ray luminosity in bins 
ofstellar mass and use the semi-analytic model of lGuo et alJ 
(l201lll to infer the underlying Lx~Ms 00 relation. Note that 
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Figure 6. The predicted 2 = 0 soft X-ray (rest-frame 0.5 —2.0 keV) luminosity—, spectroscopic temperature—, and hot gas mass-halo mass 
mean relations for the different cosmo-OWLS runs. The cyan points in the left panel represent the Lx A /500 relation oi l Anderson el all 
||2015T ) inferred from a stacking analysis of RASS data of the SDSS LBG sample. There is a clear steepening of the luminosity—halo mass 
and gas mass-halo mass relations due to gas ejection by AGN feedback. The simulations with AGN feedback reproduce the observed 
Z/X-M 500 relation well. 


LBGs are well resolved with ROSAT (i.e. #psf <C # 500 ) and, 
in any case, the X-ray luminosity is centrall y-concentrated 
due to the p 2 scaling of the X-ray emission. lAnderson et al.l 
(120151) find a significantly steeper than self-similar relation 
and conclud e that i n order to reconcile their results with 
those of the IPIntXj study, the hot gas distribution must 
become increasingly extended (less concentrated) with de¬ 
creasing halo mass. This is qualitatively consistent with the 
behaviour predicted by our simulations with AGN feedback. 

To be more quantitative, we plot in Fig.[6]the predicted 
mean soft X-ray (rest-frame 0.5 —2.0 keV) luminosity-, spec¬ 
troscopic temperature-, and hot gas mass-halo mass rela¬ 
tions at 2 = 0 for the various cosmo-OWLS models. Large 
differences are present in the predicted luminosities and hot 
gas masses at low to intermediate masses between the mod¬ 
els which do and do not include AGN feedback. 

In the left-hand panel of Fig. [S] we compare to 
the observed s t acked X-ray luminosity-mass relation of 
lAnderson et al.l (l2015lh For this comparison, we have scaled 
the observed relation to 2 = 0 assuming self-similar evolu¬ 
tion, using the mean redshift of each of the stell ar mass bi ns. 
Consi stent with the comparisons made in iLe Brun et al.l 
(12014) , we find that the observed Lx-M$oo related is well 
reproduced, both in terms of slope and zero point, by our 
AGN feedback models, with the AGN 8.0 and AGN 8.5 models 
effectively bracketing the data. 

In any case, regardless of whether the simulations match 
the observations or not, our findings highlight an impor¬ 
tant general conclusion: fluxes measured within large radii 
cannot reliably constrain the pressure profiles or integrated 
Y values within smaller radii when using matched filter¬ 
like techniques, because the obtained values depend directly 
on the assumed template profile. On the other hand, aper¬ 
ture photometry methods with aperture sizes approximately 
matched to the experiment under consideration should still 
work for these purposes (though likely at decreased signal- 
to-noise), so long as one does not try to extrapolate to 
smaller apertures. Planck has proven to be exceptionally 
good at detecting large amounts of hot gas around mas¬ 


sive galaxies (and has pushed much further down the halo 
mass scale than previous SZ experiments), thus confirming a 
fundamental prediction of cosmological structure formation. 
The next step, however, is to map the distribution of the hot 
gas in detail and for this higher resolution experiments are 
required. 

Finally, in terms of cosmological analyses with current 
Planck data, given that the flux measured within the larger 
radius 5rsoo can be reliably inferred independent of the de¬ 
tails of the assumed spatial template, and that the simu¬ 
lations all predict very similar Y 5 r500 -Msoo relations nearly 
(although not completely) independent of the sub-grid de¬ 
tails, it raises the question of whether the focus should shift 
to larger aperture measurements. If the goal is to lose sen¬ 
sitivity to the effects of uncertain baryonic physics in order 
to minimize systematic errors in cosmological analyses, this 
would indeed seem to be a sensible route to take in the fu¬ 
ture. One option would be to measure the tSZ flux within a 
large fixed physical aperture (to avoid covariances between 
the flux and mass) and to use simulations and/or a subset 
of the observations with very high quality data to calibrate 
the tSZ flux-mass relation needed for cosmological studies. 
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APPENDIX A: RESOLUTION STUDY 

We examine the sensitivity of our results to numerical res¬ 
olution. As currently available hardware prevents us from 
running higher resolution simulations in 400 h Mpc on a 
side boxes, we use smaller simulations for testing numerical 
convergence and extending the mass range of the presented 
study (see Section ^. 111 . They are 100 /i -1 Mpc on a side and 
use 2 x 256 3 particles (which is the same resolution as our 
2 x 1024 3 particles in 400 h” 1 Mpc box runs) and 2 x 512 3 
particles (i.e. eight times higher mass resolution and two 
times higher spatial resolution). They assume the WMAP7 
cosmology. Note that the convergence tests are made using 
the true physical properties of the simulated systems (i.e. no 
synthetic observations were used). 

In Fig. EH we compare the median integrated tSZ flux- 
A /500 relations at 2 = 0 at the resolution of the production 
runs (dashed lines) and at eight times higher mass resolution 
(solid lines) for the four models used here (NOCOOL, ref, 
agn 8.0 and agn 8.5). We find that the integrated tSZ flux 
is adequately converged down to M 500 ~ 2 x 10 13 M@ at the 
resolution of the cosmo-OWLS production runs, justifying 
the transition between the high-resolution simulations and 
the production runs at M 500 = 2 x 10 13 Mq adopted in the 
present study. 


APPENDIX B: SOURCE CONFUSION 


The issue of source confusion in cluster surveys, and more 
specifically in tSZ surveys (for which source confusion is 
expected to be more problematic since the strength of the 
tSZ effect is redshift independent), especially those using in¬ 
struments with rather large beams (the Planck survey falls 
into that category with angular resolutions ranging from 5’ 
to 31’ depending upon the channel frequency) has already 
been explored using simulations of various degrees of real¬ 
ism over the past fifteen years (e.g. IVojLEvrar^fc^r^arJ 


2001; IWhite. Hernquist fc Springell 120021 ; Hallman et al.l 


gel 200J; 

20071 : iHolder, McCarthy fe Babull I2007M . In brief, they all 
concluded that confusion will be an issue for tSZ surveys of 
galaxies, grou ps and low- mass clusters (M 500 < 10 14 Mq). 
For instance, IVoit et ahl (|200lh found using a back of the 
envelope calculation (using the Press-Schechter formalism) 
that the probability that any given line of sight will en¬ 
counter a virialised structure with ksT > 0.5 keV is of order 
unity and that the virialised regions of groups and clusters 
cover over a thir d of the sky. Their us e of the Hubble volume 
simulations fe.g. jjenkins et akll200lh . which has the advan¬ 
tage over the analytic calculation of taking into account the 
clustering of virialised objects, corroborated their estimate 
of the group and cluster covering factor. 
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Figure Al. Effect of numerical resolution on the median inte¬ 
grated tSZ flux— Af/500 relation at 2: = 0 . The median integrated 
tSZ flux is adequately converged down to M500 ~ 2 x 10 13 Mq, 
thus justifying the transition between the high-resolution simula¬ 
tions and the production runs at M500 = 2 x 10 13 Mq. 
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Figure Bl. The mean separation of bright back¬ 
ground/foreground sources from the central halo (i.e. the 
object of interest) as a function of the central object’s halo mass 
M 500 • The angular separation is normalized by the angular size 
0500 of the central object. The central objects are located at 
2 ^ 0.4. The solid, dashed and dot-dashed curves correspond 
to the distance to objects with Y 500 > T 50Q cen t ra i object/^’ 

T 500 > ^500,central object an( f ^500 > ^ 500 ,central object’ 
respectively. If the fluxes are estimated within 5r5oo> confusion 
is expected to become important for individual haloes with 
M500 < 10 14 Mq. 
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Figure B2. Testing the effects of source confusion on the recov¬ 
ered mean Y 500 -M 5 Q 0 relation. The blue dashed and cyan dot- 
dashed lines correspond to the e mpir ical best-fitt ing (see their 
equation 22) scaling relations of lAld and IPIntXj . respectively. 
The fill ed cyan circles with error bars represent the observational 
data of I PintXll . The filled orange circles correspond to the un¬ 
weighted mean input Y 500 ( comp uted using the empirical best¬ 
fitting Ysqq-Msqq relation of lAlOl) . The empty green squares and 
magenta triangles with error bars correspond to the weighted 
mean and errors for the maps generated by injecting the GNFW 
haloes corresponding to the input Y 500 at random positions with 
and without using the original tSZ maps from the NOCOOL simu¬ 
lation as background. The mean recovered flux is relatively insen¬ 
sitive to the presence of a tSZ background when a large number 
of haloes are stacked. 


In Fig. EH we show the mean angular separation 6 nor¬ 
malized by the angular size 6500 of the central object as a 
function of the central object’s A/500 • The central objects 
are located at z ^ 0.4 in order to mimic the LBG sam¬ 
ple. The solid, dashed and dot-dashed curves correspond 
to the distance to objects with Y500 > Y 500 cen t ra l object/^’ 

Y 500 > ^ 500 ,central object ant ^ ^ 500 > ^ 500 ,central object’ 
respectively. The integrated tSZ flux Y 500 was computed 
approximately using the self-similar prediction: Y 500 oc 
E(z) 2 / 3 AfgQQ / d\ where </a is the angular diameter distance. 
Both panels use the high-resolution version of the AGN 8.0 
simulation below A /500 = 2 x 10 13 M 0 and the production 
run above that threshold. Confusion is hence expected to 
be a source of large scatter (it could be more than a 100 
per cent) at the low-mass end for estimates of the tSZ flux 
of individual haloes, as most of them have an object with 
a similar tSZ brightness (either in the foreground or the 
background) which overlaps with them in projection (i.e. 
8/8500 < 2). Note that this can lead to both flux over¬ 
estimation and underestimation as haloes with an over¬ 
lapping neighbour with a comparable tSZ flux will have 
their flux boosted (could be up to doubled), whereas haloes 
with no overlapping neighbour with a comparable tSZ flux 


could have their flux underestimated (or even become unde¬ 
tectable) due to background overestimation. 

While confusion is expected to result into sizeable errors 
in the recovered tSZ fluxes of individual low-mass haloes, 
the recovered mean Sunyaev-Zel’dovich flux (from stacking 
a large number of systems in mass bins) can still be un¬ 
biased. To test this hypothesis, haloes that were generated 
usin g a template based on the universal pressure profile of 
were injected into the original Compton y maps. The 
results of this test are presented in Fig. IB2I For each of the 
haloes in the synthetic LBG catalogue of the NOCOOL simu¬ 
lation, the normalization of the profile (or equivalently I500) 
was set using the empirical Y500 — A/500 relation of|Aly (see 
their equation 22; dashed blue line). A new position was 
then drawn at random and the flux distributions from all 
the haloes in the catalogue were combined in order to gen¬ 
erate a new Compton y map on which the MMF was run 
both with and without using the original Compton y maps 
of the NOCOOL simulations as background. The filled orange 
circles correspond to the unweighted mean input Y500 in 
bins of A/500 ■ The empty green squares and magenta trian¬ 
gles with error bars correspond to the weighted mean and 
errors for the maps gene rated by injecting the generalized 
NFW (iNasai et al.ll2007l: see Section 0 haloes correspond¬ 
ing to the input Y500 at random positions with and without 
using the original thermal tSZ maps from the NOCOOL sim¬ 
ulation as background. In order to limit the impact of the 
added noise (see Section [2.311 . a 100 realisations of the small 
higher realisation simulation had to be used (instead of the 
usual 25) and 25 realisations (instead of the usual 2) were 
used for the large lower resolution simulation. The obtained 
results are consistent with there being no effects of uncor¬ 
related confusion over the whole mass range (as the ma¬ 
genta and green symbols are consistent within their error 
bars with their input values) when a large number of haloes 
are stacked. The mean recovered flux is thus relatively insen¬ 
sitive to the presence of a tSZ background. This insensitivity 
to uncorrelated structures could be due to the fact that the 
MMF does not use the k = 0 mode (which corresponds to 
the constant mean y value of the map). 


APPENDIX C: Y 5R 500 TO Y500 CONVERSION 
FACTOR 

In Fig. Ell we compare the conversion factor between 
Y5r500 and Y500 derived (numerically) from the new mass- 
dependent spatial template with the true mean relation de¬ 
rived directly from the simulations (i.e. using the full 3D 
particle distribution). Overall the template reproduces the 
true trend well, but tends to slightly overestimate (underes¬ 
timate) the conversion factor at the highest (lowest) masses. 
Note that these small differences are apparent when one 
compares Fig. 0 and Fig. 0 (left panel) in detail; compare 
the recovered fluxes (empty black diamonds) to the true re¬ 
lation (solid black curve) in the two plots and it is apparent 
that the inferred level of bias in the recovered flux differs 
slightly between the two plots. 

This paper has been typeset from a Tj^X/ DTj^X file prepared 
by the author. 
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Figure Cl. Comparison of the conversion factor between YsrSOO 
and y 50 o derived from the new mass-dependent spatial tem¬ 
plate (solid curve) with the true mean relation derived directly 
from the simulations (dashed curve). The horizontal dotted line 
sh ows the conversion factor for the ‘universal pressure profile’ 
of lArnaud et al] (|2010h . Overall the template reproduces the true 
trend well, but tends to slightly overestimate (underestimate) the 
conversion factor at the highest (lowest) masses. The universal 
pressure profile significantly underestimates the magnitude of the 
conversion factor for galaxies and groups. 
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